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Domain walls in superconductors: Andreev bound states and tunneling features 

S. P. Mukherjee and K. V. Samokhin 

Department of Physics, Brock University, St. Catharines, Ontario, Canada L2S 3A1 

Domain walls can be formed in superconductors with a discrete degeneracy of the ground state, for 
instance, due to the breaking of time reversal symmetry. We study all cases where the formation of 
domain walls is possible in a tetragonal superconductor with the point group D4/1. We discuss both 
triplet and mixed singlet order parameters. It is found that in all cases the domain walls support 
subgap Andreev bound states, whose energies strongly depend on the direction of semiclassical 
propagation. The bound state contribution to the density of quasiparticle states exhibits peculiar 
features, which can be observed in tunneling experiments. 

PACS numbers: 74.20.-z, 74.55.+V 


I. INTRODUCTION 

The order parameter inhomogeneities in unconven¬ 
tional superconductors and superfluids, such as the 
Abrikosov vortices, twin and grain boundaries, interfaces, 
magnetic impurities, and domain walls (DWs) are the 
breeding ground for gapless localized quasiparticles with 
distinct topological properties.^ On the other hand, the 
very presence of DWs is a direct evidence of an unconven¬ 
tional nature of the pairing, because the DWs can appear 
only in a superconductor with two or more distinct de¬ 
generate ground states, which transform one into another 
by some discrete symmetry operations, e.g. by time re¬ 
versal. This is possible if the superconducting order pa¬ 
rameter has more than one component, i.e. corresponds 
to either a multi-dimensional irreducible representation 
of the crystal point group or to a mixture of different, 
e.g., one-dimensional (ID), representations. 

One of the best-studied examples of superconductors 
or superfluids which can support DWs is a spin-triplet 
chiral p-wave state, which might be realized in Sr2Ru04 
(Ref. 2) or in thin films of superfluid ^He-A (Ref. 3). 
The order parameter in the chiral p-wave state corre¬ 
sponds to the odd two-dimensional (2D) representation, 
Eu, of the tetragonal point group D^h and is propor¬ 
tional to kx ± iky. DWs separate the opposite chirality 
states, kx + iky and kx — iky, which break time rever¬ 
sal symmetry (TRS) and have the same energy in the 
absence of external magnetic field. One possible mech¬ 
anism of the DW formation is that domains of opposite 
chirality appear spontaneously in different parts of the 
system upon cooling across the phase transition, due to 
the sample inhomogeneity. Alternatively, an increase in 
the gradient energy near the DW might be compensated 
by the creation of low-energy quasiparticles bound to the 
DW, which is particularly effective in one-dimensional 
systems."^ Experimentally, there is evidence of DWs in 
Sr2Ru04, both from the magnetic field modulation of 
the single-face Josephson critical current® and from the 
anomalous hysteresis of voltage-current characteristics,® 
and also in slabs of superffuid ®He (Ref. 7). Nonchiral 
states with the order parameters proportional to kx ± ky 
or kx , ky are also phenomenologically possible for p-wave 


pairing.® In both cases the ground state is twofold degen¬ 
erate and, therefore, DWs can exist. 

Other examples of superconducting systems which 
admit the existence of DWs are d + is and d + id 
states. In these cases the pairing is singlet and pre¬ 
dominantly d-wave, with an admixture of a subdomi¬ 
nant component of a different symmetry, s-wave or an¬ 
other d-wave, in the bulk. The possibility of such fully 
gapped states, which break both time reversal and a dis¬ 
crete lattice symmetry, had been first discussed theoret¬ 
ically for high-Tc cuprates,® but it was later ruled out 
based on phase-sensitive tests of pairing symmetry.^® Re¬ 
cently, there has been a renewal of interest in the d + is 
states in the context of iron-based superconductors, both 
theoretically^^ and experimentally.^® Various other TRS- 
breaking states have also been discussed for URu2Si2 
(Ref. 13), Bai_^K,,Fe2As2 (Ref. 14), UPta (Ref. 15), 
doped graphene (Ref. 16), undoped bilayer silicene (Ref. 
17), SrPtAs (Ref. 18), Naa;Co02 • yH20 (Ref. 19), and 
PrOs4Sbi2 (Ref. 20). 

Motivated by the growing evidence of the existence of 
multiply degenerate, in particular TRS-breaking, states 
in unconventional superconductors, in this paper we 
study how the DWs separating different ground states 
affect the spectrum of fermionic quasiparticles. The or¬ 
der parameter variation across a DW creates a quasi¬ 
particle bound state, resulting in characteristic features 
in the density of states (DOS), which can be detected, 
e.g., in tunneling measurements. We focus on the quasi- 
2D tetragonal case, which is applicable to both Sr2Ru04 
and the iron-based materials. While most of the previous 
works considered only the DW bound state spectrum for 
the chiral p-wave state,we also study the d + is and 
d + id states, as well as nonchiral p-wave states. While 
we consider mixed pairing states in the bulk, another 
possibility is that a subdominant order parameter is sta¬ 
bilized only near a surface, in which case the spectrum 
of the surface bound states was found in Ref. 23. The 
paper is organized as follows. In Sec. II, we derive a 
general expression for the DW bound state energy in the 
semiclassical (Andreev) approximation. In Sec. Ill, we 
list all possible cases which admit the DW formation for 
the point group D^h- In Secs. IV and V, we discuss the 
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FIG. 1 : (Color online) Schematic diagram showing the DW, 
the semiclassical trajectory of a Bogoliubov quasiparticle, and 
the order parameters A± on both sides. 

triplet and mixed singlet cases, respectively, and calcu¬ 
late the spectrum of the bound states as well as their 
contribution to the quasiparticle DOS. Finally, we sum¬ 
marize our results in Sec. VI. Throughout the paper we 
use the units in which h = e = c = 1. 


II. ANDREEV BOUND STATES 

We consider a quasi-2D superconductor with a cylin¬ 
drical Fermi surface and assume that the DW is along 
the y-axis, so that the superconducting domains are ex¬ 
tended on both sides of the DW to infinity along the x- 
axis. The order parameter variation takes place within a 
region of width near the DW. Throughout our deriva¬ 
tion we use the Andreev approximation,in which the 
Bogoliubov quasiparticles propagate along semiclassical 
trajectories characterized by the Fermi-surface wavevec- 
tors kf = {kF,x,kp^y) = kp^cosOjSmO). This approxi¬ 
mation is justified since the scales of variation of the su¬ 
perconducting order parameter, including the DW width 
^d, are much greater than the Fermi wavelength kp^. 

For a given direction of semiclassical propagation, the 
gap function takes the form Akjr(x) = A{9,x), with the 
0-dependence determined by the pairing symmetry, see 
Sec. Ill below. The unavailability of an exact analytical 
expression for the a:-dependence of the DW order param¬ 
eter leads to various approximation schemes, see Refs. 26 
and 27 for the chiral p-wave case. We use a simple, an¬ 
alytically treatable, model in which the DW is assumed 
to be a sharp boundary of zero width located at a; = 0. 
Then the gap function has the following form: 

Akp(x) = A+(0), x>0, 

Akp{x) = A_(0), a: < 0, ^ ^ 

i.e. it is characterized by two different complex numbers 
along each semiclassical trajectory, see Fig. 1. 

In the Andreev approximation, the quasiparticle wave 
function is a product of a rapidly oscillating plane wave 


gikpr ^ slowly varying envelope function 'F(x), which 
satisfies the equation 

where vf,x = vf cosO and vf is the Fermi velocity. Us¬ 
ing the sharp DW model, Eq. (1), we find that the 
bound-state wave function has the form ll'±(a;) ^ 
at a; > 0 (a; < 0), where k± = ^±/\vf,x\, with rt± = 
a/|A±P — >0. In the general pairing case, the en¬ 

velope wave function has four components, but in all the 
models we consider below the spin channels decouple and 
the equations can be reduced to a two-component form. 

Substituting the wave functions into the boundary con¬ 
dition 0) = 'I'(-I-O), we arrive at the following equa¬ 
tion for the bound state energy: 


E — iVF,xK+ E + iVF,xK-^ 
which can be represented in the form 

E+ in_ =-f{E -in+), (4) 

where E = E sgn{vF,x) and 7 = A_/A+ = jr + 17 /. 
Taking the real and imaginary parts of Eq. (4), we ob¬ 
tain: 

Qp = = hLjlJAE. ( 5 ) 

7/ 7/ 

The Hrst of these equations gives 

E = sgn(A)-=i=|A+|, (6) 

VT+W 

where 

/3 = (7) 

li 

The sign function that appears on the right-hand side of 
Eq. ( 6 ) can be assigned using the condition that D+ > 
0, from which we find sgn(E) = sgn(/3). Therefore, we 
finally obtain for the Andreev bound state (ABS) energy: 

EbiO) = \A+{0)\ ^ =sgn[/3(0)cos0]. ( 8 ) 

y/l+ 

There is an additional constraint coming from the second 
of the expressions (5). Since r2_ > 0, we have 

sgn(| 7 p- 7 fl)sgn(l- 7 fl) = 1. (9) 

Thus, for each direction of semiclassical propagation sat¬ 
isfying the condition (9), there is one ABS with the en¬ 
ergy given by Eq. ( 8 ). Eor singlet pairing, the gap func¬ 
tion is even in momentum and A±(0 -|- tt) = Aj-(0), 
while for triplet pairing, the gap function is odd and 
Aj-(0 -I- tt) = —A±{6). It follows from Eqs. (7) and ( 8 ) 
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that /3(0 + 7 r) = j3{9) and, therefore, Eb{9 + -K) = —Ei,{9) 
in all cases. 

The ABS wave function is exponentially localized on 
both sides of the DW, with the characteristic scales given 
by k^^{9) = |n_F,x|/\/|A±p — E'^. The sharp DW ap¬ 
proximation is valid for those directions of semiclassical 
propagation for which the DW width is smaller than 
This condition is strongly angle-dependent and, in 
particular, fails for “grazing” trajectories, for which 9 
is close to ± 7 r/ 2 , therefore, vf,x 0 and the Andreev 
approximation itself is inapplicable. However, for most 
directions of fep, one can use the following rough esti¬ 
mate: > uf/Aq where Aq is a characteristic 

value of the gap and ^ is the superconducting correlation 
length. Therefore, the sharp DW model is quantitatively 
justified if ^ ^ ^±^- 

The ABS’s contribute to the quasiparticle DOS, which 
can be probed in tunneling experiments. The momen¬ 
tum along the y-axis is conserved and the bound-state 
DOS per unit length of the DW is given by the following 
expression: 

N,{E)= r '^Y.8[E - Ed{\^F)l 

where kp^x = A “ ^Fy Making the change of vari¬ 
ables kp^x = kp COS0, kp^y = kp sin0, we obtain: 

p2-K 

Nb{E)=Npvp d9 \ cos9\6[E - Eb{9)], (10) 

Jo 

where IVf is the Fermi-surface DOS in two dimensions. 


III. NONUNIFORM STATES IN A 
TETRAGONAL SUPERCONDUCTOR 

The point group has ten irreducible representa¬ 
tions in total (considering even as well as odd representa¬ 
tions), of which eight are ID and two are 2D. The forma¬ 
tion of DWs is possible only for those superconducting 
classes which are degenerate with respect to some dis¬ 
crete symmetry.®’^® We consider the cases of singlet and 
triplet pairing separately. 

Triplet order parameter is described by a spin vector 
d(k). We assume that the spin-orbit coupling fixes the 
direction of d(k) along the z axis. Since the ID rep¬ 
resentations cannot support DWs, we consider the odd 
irreducible representation E^, which corresponds to p- 
wave pairing. The semiclassical gap function takes the 
following form: 

AkF(a:) = ili{x)kp^x + V2{x)kp,y, (11) 

where r] = (??i,? 72 ) is the two-component order param¬ 
eter, which transforms like a vector in the xy plane, 
and the basis functions of E^ are given by kp^x = cos 9, 
kp^y = sin0. In Table 1 we show all stable states for 


TABLE 1: Superconducting states corresponding to the irre¬ 
ducible representation Eu of the point group D^h with strong 
spin-orbit conpling. The hrst column lists the possible su- 
percondncting classes, where R stands for time reversal (the 
notations are the same as in Ref. 26), the second and third 
columns show the number of degenerate ground states and 
the specific forms of the order parameter in the two domains, 
and the vector d(k) is given in the last column. 


Dih 

degeneracy 

■n 

d(k) 

DiiE) 

2 

(lA) 

z(kx -1- iky) 



(1,-*) 

iky ) 

Di{C2) X R 

2 

(1,0) 

zkx 



(0,1) 

zky 

D2{a)xR 

2 

(1,1) 

^{kx “1“ 



(1,-1) 

z{kx - ky) 


which the formation of DWs is possible. Note that only 
the states with 77 oc (1,±*) break TRS and do not have 
gap nodes in the bulk, while the other two states pre¬ 
serve TRS but break some discrete crystal symmetries 
and have bulk gap nodes. 

In the singlet case and quasi-2D geometry (the latter 
implies fcz-independent basis functions), TRS breaking 
and DWs are possible only in the states formed by a 
combination of two pure states belonging to two different 
ID irreducible representations.^® The semiclassical gap 
function can be written as 

= ll)i{x)Lpi{\cp) + 'tl>2[x)Lp2(}i-F), ( 12 ) 

where il^i ^2 are the order parameters and ipi ^2 are the 
basis functions of the ID representations, satisfying 
<Pi, 2 (kF) = <Pi, 2 (—kF). We consider the following three 
possibilities: a mixture of Big and i? 2 g representations, 
for which 

(pi — kpj. kpy-f ip2 — ‘^kpxkp^y ( 1 ^) 

a mixture of Big and Aig representations, for which 

‘Pi =kF^x-kly, = (14) 

and a mixture of i? 2 g and Aig representations, for which 

ipi = 2kF^xkF,y, P2 = 1. (15) 

The superconducting state in the bulk is given by 
('i/'i,'!/) 2 ) oc (Ai,±iA 2 ), where Ai and A 2 characterize 
the strengths of the individual pairing components. The 
first of the above states is referred to as a dx 2 -y 2 ± idxy 
state, while the last two are dx^-y^ ± is and dxy ± is 
states, respectively. All of these states are fully gapped 
in the bulk.^® 


IV. TRIPLET STATES 

In this section we discuss the ABS spectra and the 
corresponding DOS for the three triplet p-wave states 
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FIG. 2: (Color online) The variation of the ABS energy with 
the direction of the quasiparticle propagation for domains 
with rj oc (1, ±i), for different values of X- 


listed in Table I. 



FIG. 3: (Color online) The energy dependence of the ABS 
contribution to the DOS (in arbitrary units) for domains with 
r; oc (1, ±i), for different values of X- 


A. T 7 oc (1, ±i) 


all 0 and x- For the ABS energy we then obtain: 


In the chiral p-wave state, the order parameters on 
both sides of the DW have the following form: 


Ebje) 

Aq 



cos 0 sin 



(17) 


?7 = Ao(l,i), x<0, 

T] = Ao(l, -i)e*^, a: > 0, 

where Aq is the gap magnitude and x is the common 
(Josephson) phase difference between the domains. The 
latter has to be included in order to satisfy the current 
conservation across the DW. Its value depends on the mi¬ 
croscopic details of the system, see Appendix A. Without 
loss of generality we assume 0 < x < tt. 

It follows from Eqs. (1), (7), and (11) that A_(0) = 
Aoe*®, A+(6») = Aoe*Xe-*®, and 

13 = tan (^ - I) ■ 

Note that the condition (9) is automatically satisfied, 
since \x\^ — XR = 1 “ a-nd Eq. (8) can be used for 


see also Ref. 22. In Fig. 2 we plot this last expression for 
three representative values of the common phase differ¬ 
ence X- The ABS energy has discontinuities at the spe¬ 
cial directions of semiclassical propagation: at 0 = ±7r/2, 
i.e. for the quasiparticles moving parallel to the DW, in 
which case the Andreev approximation is not applica¬ 
ble, and also at 0 = x/2 and 0 = x/2 -b tt, for which 
A+ = A_ and, therefore, the DW is “invisible” to the 
quasiparticles. 

Next, we calculate numerically the ABS contribution 
to the DOS per unit length of the DW by using Eq. (10) 
and plot it in Fig. 3. One notable feature of them is the 
presence of van Hove singularities at the gap edge, i.e. 
at Eb = ±Ao for x tt. These singularities are due to 
the vanishing of the slope of Eb(0) at 0 = x/2 and 0 = 
x/2 + TT. For X = ^ tFe slope vanishes for 0 = tt/2, 37r/2, 
where cos 0 = 0 and the van Hove singularities are absent. 
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FIG. 4: (Color online) The variation of the ABS energy with 
the direction of the quasiparticle propagation for domains 
with Tj oc (1,0) and (0,1), for different values of x- Regions 
bounded by blue vertical dotted lines and not containing any 
bound states correspond to the directions of propagation for 
which the condition (9) is not satisfied. 

In fact, in the latter case, the DOS is independent of 
energy. 

B. r] oc (1, 0) or (0,1) 

The order parameters in this nonchiral p-wave state 
have the form 

T7 = Ao(1,0), X < 0, 

77 = Ao(0,l)e*>^, x>0, ^ 

therefore, A_(0) = Aocos0, A+(0) = Aoe*^sin0, and 

^ cos Y — tan 9 
P = -- 

smx 

We calculate numerically the ABS spectra using Eq. ( 8 ) 
for different representative values of x present the 
results in Fig. 4. Note that for some directions of semi- 
classical propagation the condition (9) is not satisfied 
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FIG. 5: (Color online) The energy dependence of the ABS 
contribution to the DOS (in arbitrary units) for domains with 
rj oc (1, 0) and (0,1), for different values of y. For y = 0 the 
DOS has a delta-function singularity at i? = 0. 


and, therefore, there are no bound states. The case of 
y = 7 r/ 2 , when 7 ^; = 0, is an exception: the ABS’s exist 
at all 9. 

The ABS contribution to the DOS is plotted in Fig. 5. 
For y = 0, the DOS has a delta-function singularity at 
zero energy due to the flatness of the corresponding ABS 
dispersion. This can be understood as follows. At y = 
0, the gap functions on both sides of the DW are real, 
A_(0) = AqCOS^ and A+(0) = Apsinfl. If sin 20 < 0 
then A+ and A_ have opposite signs. In this case, the 
Andreev equations are equivalent to the Witten model 
of supersymmetric quantum mechanics and the ABS can 
be shown to have exactly zero energy, see Ref. 25 and 
the references therein and also Ref. 30. The condition 
sgn(A 4 .A_) = —I is satisfied only at 7 r /2 < 9 < tt and at 
3tt/2 < 9 < 2tt, where the ABS energy is zero, as shown 
in the top panel of Fig. 4. 

Adding the phase y makes the gap functions A± in¬ 
trinsically complex, in which case the analogy with the 
supersymmetric quantum mechanics breaks down and a 
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FIG. 6: (Color online) The variation of the ABS energy with 
the direction of quasiparticle propagation for domains with 
r; oc (1,±1), for different values of X- Regions bounded by 
blue vertical dotted lines and not containing any bound states 
correspond to the directions of propagation for which the con¬ 
dition (9) is not satisfied. 


simple physical explanation of the spectra is no longer 
available. Our results show that as x increases, the ABS 
bands acquire curvature and also the range of the ABS 
existence gets broader until, eventually, the internal gaps 
completely disappear for x = '^/2, see the middle and 
bottom panels of Fig. 4. At 0 < x < 7r/2, there are two 
sets of the van Hove singularities in the DOS, correspond¬ 
ing to the vanishing of the dispersion slope at the ABS 
band centers, and also two internal gaps, corresponding 
to the absence of the ABS at certain energies. Finally, at 
X = 7r/2, the four van Hove singularities merge into two. 

C. r] (X (1, ±1) 

The order parameters in this nonchiral p-wave state 
have the form 

r 7 = Ao(l,l), a;<0, 

T] = Ao(l, a; > 0, 


- 

X=0 
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- 



FIG. 7: (Color online) The energy dependence of the ABS 
contribution to the DOS (in arbitrary units) for domains with 
r; oc (1, ±1), for different values of x- For x = 0 the DOS has 
a delta-function singularity at A = 0. 


therefore A^(0) = Ao(cos0 -I- sin0), A+(0) = 

Aoe®^(cos0 — sin0), and 

^ cosx — tan(7r/4 — 0) 
sinx 

The ABS spectra and contribution to the DOS are cal¬ 
culated numerically for different representative values of 
X and plotted in Figs. 6 and 7, respectively. 

Similarly to the previous case, the condition (9) is not 
satisfied for some values of 9 (except at x = 7'‘/2), re¬ 
sulting in the spectra consisting of several disconnected 
parts and the gaps in the DOS. The van Hove singu¬ 
larities are again present, degenerating for x = 0 into 
a delta-function singularity at zero energy. The latter’s 
origin is explained by the fact that at x = 0 the gap 
functions A+ and A_ are real and have opposite signs 
if cos 29 < 0, in which case the ABS has zero energy for 
7r/40 < 37r/4 and 57r/40 < 77r/4, as shown in the top 
panel of Fig. 6. At nonzero Xi and A_ become com¬ 
plex, which results in the ABS bands acquiring dispersion 
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FIG. 8: (Color online) The variation of the ABS energy with 
the direction of quasiparticle propagation for the domains 
d^ 2 _y 2 izidxy, for different values of x and p. Thin red line is 
for p — 0.2 and thick blue line is for p = 0.8. For x = 0 l-he 
curves superpose on each other. 


and spreading to other angles of propagation until filling 
the entire angular range at x = 



E/A 


max 


FIG. 9: (Color online) The energy dependence of the ABS 
contribution to the DOS (in arbitrary units) for the domains 
d^ 2 _y 2 ± idxy, for different values of x and p. Thin red line is 
for p = 0.2 and thick blue line is for p = 0.8. For x = 0 they 
superpose on each other. 


A., dj.2 _ y2 di idjjy 


V. MIXED SINGLET STATES 

In this section we discuss the ABS spectra for TRS- 
breaking mixtures of different singlet states. In all three 
cases, see Sec. Ill, the order parameters on both sides of 
a sharp DW can be written as 

('0i,V’2) = (Ai,iA2), x < 0, .20) 

(V'1,'02) = (Ai,-fA 2 )e*^, a; > 0, 

where Ai,A 2 > 0 are constant amplitudes of the two 
singlet channels (without loss of generality, we assume 
that A 2 < Ai) and x is the common phase difference, 
see Appendix B. 


In this case, we obtain from Eqs. (I), (12), and (13): 

A_(0) = Ai cos20 + iA 2 sin20, , . 

A+(0) = (Ai cos 20 — iA 2 sin20) 

We see that the gap magnitude is the same in both do¬ 
mains and given by 

|A+| = |A_| = Ab„/fe(0) = Af cos^ 20 -b A| sin^ 20, 

which attains its maximum value, A^ax = Ai, at 0 = 
m:l2, and its minimum value, Amin = A 2 , at 0 = (2n -|- 
l)7r/4, where n = 0,1, 2,3. Thus, the bulk gap has an 
anisotropic magnitude without nodes, while the phases of 
the gap function are different in the two domains: A_ = 
Abuike‘* and A+ = Abnike‘^e~'^'>’, where 

(piO) = arctan(/9tan20) 


























































and 


P=^, 0<p<l. (22) 

From Eq. (7) we obtain: 

/? = tan - I) , (23) 


therefore, the expression (8) for the ABS spectrum takes 
the following form: 


Eb{0) 


= + (1 — p^) cos^ 29 cos {^(j) ~ 


X sgn 


cos0sin^0—, (24) 


which is shown in Fig. 8 for different representative val¬ 
ues of p and X- Although the ABS exist for all directions 
of semiclassical propagation, the spectra are not contin¬ 
uous, due to the presence of the sign functions in Eq. 
(24). In some cases, relatively simple explicit expressions 
for the ABS energy can be obtained: for y = 0, we have 

^ = cos20sgn(sind), 

^max 


which does not depend on p, and for y = tt, 

= —psin20sgn(cos0cos20). 

^max 

The DOS, shown in Fig. 9, has the van Hove singularities, 
corresponding to the vanishing of the dispersion slope. 
For X = 0, the singularities occur at \E\ = 1, while for 
X = TT - at |f;| = p. 



FIG. 10: (Color online) The variation of the ABS energy with 
the direction of quasiparticle propagation for the domains 
d^ 2 _y 2 ± is, for different values of x and p. Thin red line 
is for p = 0.2 and thick blue line is for p = 0.8. 


and p is given by Eq. (22). Therefore, we obtain for the 
ABS spectrum: 


B. d„, 2 _y 2 ± is 

According to Eqs. (1), (12), and (14), the supercon¬ 
ducting domains are described by the following gap func¬ 
tions: 


Eb{e) 

1 + cos^ 20 j 

'6 

Amax 

V p" 

-b 1 ' 

. 2/ 


X sgn 

cos 0 sin {^(f) 

-f)] 


(26) 


A_(0) = Ai cos 20 -I- iA 2 , 

A+(0) = (Ai cos20 — iA 2 ) e*^. 

The gap magnitudes are the same in both domains: 
|A+| = |A_| = AbuikiO) = ^Al cos2 20 + Ai. 


which is plotted in Fig. 10. One can see that the ABS 
exist for all 0, but their spectra are not continuous. The 
DOS is shown in Fig. 11. In particular, in the case of x = 
TT, we find from Eq. (26) that \Eb\IAmax = p/\/p^ -I- 1, 
i.e. does not depend on 0, which results in the delta- 
function singularities in the DOS. 


The bulk gap magnitude is anisotropic, without nodes. 
It attains its maximum value, Amax = \/A^ -I- A|, at 
0 = n7r/2, and the minimum value, A^^ = A 2 , at 0 = 
(2n -I- l)7r/4, where n = 0,1, 2, 3. 

The parameter /3 is given by Eq. (23), where 


(j) = arctan 


V cos 20 


) 


C« dxy I tS 

According to Eqs. (1), (12), and (15), the supercon¬ 
ducting domains are described by the following gap func- 


















































9 



- 




- 

- 




— p=0.2 

- 

- 




— p=0.8 

- 

- 




- 


- 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 


E/A 


max 


FIG. 11: (Color online) The energy dependence of the ABS 
contribution to the DOS (in arbitrary units) for the domains 
d^ 2 _y 2 ± is, for different values of x a-nd p. Thin red line is 
for p = 0.2 and thick blue line is for p — 0.8. For x = the 
DOS has delta-function singularities. 


tions: 


A_ (0) = Ai sin 29 + 
A-|_(0) = (Ai sin 20 — il^ 2 ) 


(27) 


The gap magnitudes are the same in both domains: 

|A+| = |A_| = Afc„zfe(0) = ^Afsin^ 20 tA^. 

The bulk gap magnitude is anisotropic, without nodes. 
It attains its maximum value, A^ax = \/A^ -|- A|, at 
0 = (2n + l)7r/4, where n = 0,1, 2,3, and the minimum 
value, Amin = A 2 , at 0 = wk 12, where n = 0,1, 2, 3. 

The parameter /3 is given by Eq. (23), where 


(\) = arctan 



and p is given by Eq. (22). The ABS spectrum has the 



FIG. 12: (Golor online) The variation of the ABS energy 
with the direction of quasiparticle propagation for the do¬ 
mains dxy ± is, for different values of x and p. Thin red line 
is for p — 0.2 and thick blue line is for p — 0.8. 


following form: 


Eb{e) 


I + sin^ 20 

p2 + l 


CO8(0-|) 


X sgn 


cos 0 sin ((h — 


2/J ’ 


(28) 


which is shown in Fig. 12. Again, the ABS exist for 
all 0 and their spectra are not continuous. The DOS is 
shown in Fig. 13. In particular, in the case of x = 
we find from Eq. (28) that \Eb\/Amax = p!\/9^ + 1, ke. 
does not depend on 0, which results in the delta-function 
singularities in the DOS. 


VI. CONCLUSION 

We have shown that the DWs separating degener¬ 
ate ground states in unconventional superconductors can 
trap fermionic quasiparticles and create Andreev bound 
states with energies inside the bulk gap. Using as an 
example a quasi-2D tetragonal superconductor, we have 
studied both chiral and nonchiral p-wave states, as well 
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ated delta-function peaks in the DOS is explained by the 
sign change of the real gap function along the semiclassi- 
cal trajectory. In the mixed singlet cases, the spectra of 
d-\-is states can have four delta-function peaks at nonzero 
energies due to the Andreev bands with flat dispersion, 
which never happens in the d id state. We hope that 
the abovementioned features can be detected in tunneling 
experiments, thus providing a direct evidence of the DW 
presence and shedding light on the pairing symmetry of 
unconventional superconductors. 


FIG. 13: (Color online) The energy dependence of the ABS 
contribution to the DOS (in arbitrary units) for the domains 
dxy ± is, for different values of y and p. Thin red line is for 
p = 0.2 and thick blue line is for p = 0.8. For y = tt the DOS 
has delta-function singularities. 


as the mixed singlet states oi d + id and d + is symme¬ 
tries. Time reversal symmetry is broken in each case, 
except the nonchiral p-wave states. 

We derived a general expression for the ABS energy in 
the semiclassical approximation, for a sharp DW model. 
The bound states exist only for certain directions of the 
quasiparticle propagation, and there is only one bound 
state for each “allowed” direction. The ABS spectrum 
is highly anisotropic, forming the ABS “bands” as the 
direction of the Fermi wavevector varies, which results in 
some prominent features in the quasiparticle DOS, such 
as internal gaps, abrupt steps, and the van Hove singu¬ 
larities. 

Some of the DOS features are specific to certain pair¬ 
ing symmetries and can help to distinguish them. For 
instance, the ABS spectrum in the chiral p-wave state 
has neither internal gaps nor delta-function peaks, in 
contrast to the nonchiral ones. In the latter cases, the 
origin of the flat zero-energy ABS bands and the associ- 
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Appendix A: Calculation of x for the triplet states 


Superconducting current can be obtained from the gra¬ 
dient terms in the Ginzburg-Landau free energy density. 
For a p-wave triplet order parameter in a tetragonal crys¬ 
tal the latter have the following form:® 

Fgrad = + A 2 (V,r 70 *(V,P,) 

+A3(V,?7,)*(V,770 + A4(V,p,)*(V.77*), (A1) 

where the summation over repeated indices, i, j = x, y, is 
implied. Replacing the gradients by the covariant deriva¬ 
tives, V —V -f 2iA, and varying with respect to the 
vector potential A, we obtain: 

j, = 2Im(Ai77;V,p, + 

+ K4r/*V^7j,). 

The order parameter components can be written in the 
following general form: 

771 = Ao/i(x)e*'^("\ 772 = (A2) 


where /i _2 are dimensionless amplitudes, (/> is the com¬ 
mon (or Josephson) phase, and 7 is the relative phase. 
From these expressions we calculate the components of 
the supercurrent as follows: 

j, = 2Al (A 1234/1 + Kifl) V,<^-2A2Ai/|V, 7, (A3) 
where K 1234 = Ki K 2 + K 3 K 4 , and 

jy = 2A2[sin7 (A2/2 V,/i - Ag/iV^/s) 

+/ 1/2 cos 7 [(A 2 + A 3 )V,((. - A 3 V, 7 ]. (A4) 


From the current conservation we get jx = const, where 
the constant is determined by the boundary conditions 
and can be set to zero. This gives us a relation between 
the phase gradients: 


Va;</> 


Kifi 

Kl23Afi + ATi/I 


Va;7. 
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Therefore, 


X = (j){+oo) - (j){-oo) 

K,fi 


dx 


^1234/1 + Kifi 


'^xl- 


the current across the DW and along the DW, with the 
following results: 


(A5) 


jx = 2lm[Ki{V + K 2 iV oo’4’2)i’2 

+K3{ri^x^2+r2^x^i)] (B2) 


We see that the parameter x depends on the microscopic 
characteristics of the system, through the coefficients Ki, 
as well as on the DW structure, described by fi, 2 {x) and 
7(0:), which are model-specific. Below we calculate x for 
different triplet states in a simple approximation in which 
the order amplitudes are constant, either globally, i.e. at 
all X, or separately at a: > 0 and a: < 0. 

In the chiral case with r] oc (l,±i), we assume that 
/i = /2 = 1 at all X, see Ref. 26. Since 7(±oo) = ±7r/2, 
we can calculate the integral on the right-hand side of 
Eq. (A5), with the following result: 


Ki 

Ari234 + Kl 


(A6) 


For the nonchiral state with rj oc (1,±1) we have fi = 
f 2 = 1, but 7(—00) = 0 and 7(-|-oo) = tt. From Eq. 
(A5), we obtain the same expression (A6) for y. For the 
nonchiral state with rj oc (1,0) or (0,1), we have /i = 
G{—x) and /2 = 0(a;), where 0(x) is the Heaviside step 
function. In this case, the current conservation is satisfied 
if both and 7 take arbitrary constant values throughout 
the system, which means that the order parameter can 
be written as = AoO(—x), r ]2 = Ao0(a:)e*^, where 
the value of y is arbitrary in our approximation. Using 
Eq. (A4) one can show that the net current along the 
DW is nonzero in the chiral state, but vanishes in both 
nonchiral states. 


Appendix B: Calculation of y for the singlet states 

As an example, we consider only the dx 2 -y 2 ±is state, 
the calculation for dxyAis and dx 2 -y 2 ±idxy states being 
similar. The gradient part of the Ginzburg-Landau free 
energy density has the form^^ 


and jy = 0. 

Representing the order parameters in the form 

= Ai(x)e *'^("), ^2 = A2 , (B3) 

where Ai_2 are the amplitudes, </> is the common phase, 
and 7 is the relative phase, we find from the current 
conservation that the phase gradients are related: 


= 


K2AI + a:3AiA2 C 0 S 7 


KiAl 


■ K 2 AI + 2 K 3 A 1 A 2 cos 7 


V,7. (B4) 


Therefore, 


y = 0(-foo) - 0(-oo) 

K 2 AI + K 3 A 1 A 2 COS 7 


dx 


KiAl + K 2 AI + 2 K 3 A 1 A 2 cos 7 


Va;7. 


(B5) 


The relative phase satisfies the following boundary con¬ 
ditions: 7(±oo) = ±7r/2, 

Assuming constant amplitudes Ai and A2, we have 



Cl -I- cos 7 
C2 -I- COS 7 


dj, 


(B6) 


where 


K 2 A 2 KiAl -f K 2 AI 

^^“ifaAi’ 2X3A1A2 ■ 

The integral in Eq. (B6) can be solved analytically and 
we finally obtain: 



2(ci - C2) 
-1 


arctan 




C2 - 1 
C2-f 1’ 


Fgrad = Ki\V^,\^+K2\Vi;2\^ 

+A3[(VyV'i)*(Vj;V'2) - (V,V'i)*(V.V'2) + C.C.] (Bl) 

(for the d + id state the K 3 term is absent). We assume 
that the order parameters '>pi ,2 depend only on x. Follow¬ 
ing the same procedure as in Appendix A, we calculate 


for C2 > 1, and 

^ ^ (ci - C2) / C2 -I-V'cj-l ^ 

2 a/c^ - 1 \ C2 - 1 -f Vc| - 1 j 

for C2 < 1. 


















